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By Percy J. Bobbitt and Peter J. Maxie, Jr. 


SUMMARY 


In order to calculate the induced loading on a store , missile, or 
pylon situated in close proximity to a wing-fuselage combination, a 
detailed knowledge of the flow field is required. The present paper 
provides some of this needed information by presenting equations and 
charts that enable the determination of sidewash in the vicinity of semi- 
infinite triangular wings at small angles of attack. These results may 
also be used directly in sidewash determinations for the conical part of 
the flow exterior to wings having finite spans and chords. At points in 
the flow field affected by the finite-wing wake or tip, additional con- 
siderations are necessary to determine the sidewash. Both the subsonic- 
wing and super sonic -wing leading-edge conditions have been treated; hence, 
sidewash is obtainable for all supersonic Mach numbers and leading-edge 
sweep angles for sweptback wings. 


INTRODUCTION 


The ever increasing use of external stores and missiles on aircraft 
has emphasized the need for a method of predicting the loads acting on 
these stores and missiles and also on the pylons by which they are often 
attached to the wing. Store, missile, and pylon loadings are required 
in designing the pylon structure, in predicting the lateral stability 
of aircraft, in determining the jettison characteristics of stores, and 
in computing the trajectories of missiles. 

A number of experimental investigations have been conducted in the 
past several years in an attempt to gain some insight into the origin 
of store loads and, in a few isolated cases, pylon loads. Through system- 
atic tests on the effects of store position, store size, wing plan form, 
pylons of various types, and Mach number, the understanding of the various 
interference effects has been increased. However, due to the many vari- 
ables involved in the airplane- store or airplane -missile problem, there 
remains a very definite need for an analytical or semiempirical method 
capable of indicating trends (if not magnitudes) when the many variables 
involved are changed. 
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One of the most important prerequisities to making a calculation 
of the store, missile, or pylon loading is a detailed knowledge of the 
flow field in the immediate vicinity of the wing. In reference 1 an 
extensive theoretical treatment of the pertinent flow fields at subsonic 
speeds has been presented, together with some experimental verification. 

At supersonic speeds the picture is not quite so complete. The longi- 
tudinal and vertical velocity components in the vicinity of the wing 
which are due to the wing angle of attack have been derived in refer- 
ences 2 , 3 , by and and some information related to the velocity com- 
ponents in fuselage flow fields and wing-thickness flow fields may be 
extracted from references 6 and 7* The lateral flow velocity in the 
vicinity of the wing due to the wing angle of attack has not yet been 
treated. 

The present paper is concerned with the latter, that is, the analyt- 
ical determination of the lateral flow velocity (sidewash) in the vicinity 
of the wing due to the wing angle of attack. The methods used parallel 
the linearized lifting- surface methods employed in reference 3 to deter- 
mine the downwash in the flow field exterior to flat lifting triangles 
of infinite chord. The analyses apply for all leading-edge sweep angles 
for sweptback wings and supersonic Mach numbers, both the subsonic-wing 
leading-edge and super sonic -wing leading-edge conditions being considered. 
It should be pointed out that the infinite -chord triangular-wing sidewash, 
which is derived herein, may also be used directly in sidewash determina- 
tions for some regions near wings having finite spans and chords. For 
points affected by the wing wake or wing tip, additional considerations 
are, of course, necessary (see ref. 5 ) and ‘the expressions for the side- 
wash velocity given in the present paper represent only a part of the 
total sidewash. The results of the investigation are presented in chart 
form for various combinations of leading-edge sweep and Mach number. 


SYMBOLS 


The system of axes used in the analysis and the positive directions 
of the velocities are shown in figure 1. 

b span of triangular wing 

C constant of integration 

c root chord of triangular wing 

E i imaginary part of incomplete elliptic integral of second 

kind 

E r real P art of incomplete elliptic integral of second kind 
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E [\j^ ~ P 2 m 2 

E(a),k) 


F r 

F(o>,k) 


I.P. 


k 

k' 


m 


complete elliptic integral of second kind with modulus 

\fl - pSm 2 

incomplete elliptic integral of the second kind with 

oO) \jl - k^o^ 2 da^ 


modulus k and sine amplitude cd, 


0 


iF 


“L 


imaginary part of incomplete elliptic integral of first 
kind 

real part of incomplete elliptic integral of first kind 


incomplete elliptic integral of first kind with modulus k 

n CD 

and sine amplitude od, / 

Jo 


do>^ 


\ j 1 " / 


1 - o>, 2 \/l - k2a>^2 


imaginary part 


modulus of elliptic integrals, \J 1 - s 0 ^ 
comodulus of elliptic integrals, k' = ^/l - 
free-stream Mach number. 


k 2 = s Q 2 


Vo, 


Velocity of sound in free stream 


cotangent of sweep angle of wing leading edge 
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So - 

U,v,w 

V 

Voo 

V 

x,y,z 


1 - V 1 - p 2 !!! 2 


pm 


a 

p 

5 




Moo - 1 


A = sn 2 (F r ,k) 


perturbation velocities in x-, y-, and z-directions 

complex side wash function, v + iv 
free-stream velocity 
harmonic conjugate of v 

Cartesian coordinates of field point with origin at wing 
apex (the positive directions are indicated in fig. 1) 

angle of attack, radians 

imaginary part of a complex variable 
complex variable, £ = p + iq. = 


p I + ^ S 


1 + \/i - (y/x) - (z/x)‘ 


complex incomplete elliptic f met ion of first kind (real 
part is given by p r and imaginary part by up) 




complex variables 


a = sn 2 (Fj_,k ! ) 

T 


0 


real part of a complex variable 
perturbation-velocity potential 


Subscripts : 


1 

2 

i 

r 


refers to complex variable £ 

refers to complex variable t) 

refers to imaginary part of a complex function 

refers to real part of a complex function 
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PRESENTATION OF SIDEWASH EQUATIONS 


The methods used in reference 3 for the determination of the down- 
wash in the flow fields exterior to flat lifting triangular wings of 
infinite chord may also be used for the determination of the sidewash 
in the flow fields exterior to these same wings. For this reason, this 
section will be restricted primarily to the presentation of the sidewash 
equations with the actual derivations given in appendixes A and B. 

Sidewash equations presented will apply directly for points below 
the right-hand panel of the wing; however, since the sidewash is anti- 
symmetric with z (i.e., about the plane of the wing) and y, the signs 
on the equation may be changed to give the sidewash for any point in 
the field. 


The sidewash at a point (x,y, z) below the right-hand panel of the 
lifting triangular wing, as developed in appendix A, is given by 


Subsonic-Leading-Edge Wing 



/ 


H 1 K " ' ' ' 




(1) 


where 


1 -\jl - |B 2 m2 


pm 


6 


NACA TN 3938 


In order to calculate the A and a functions in equation (l), 
proceed as follows. First, calculate p and q: 



and then calculate r^, 5-^, r 2 , and 5 2 : 


T l = 


p|/ \l[v 2 - q 2 + s 0 2 y + 4p 2 q 2 + p 2 - q 2 + s Q 2 + 


p 2 - q2 + s 0 2 ' + 4p2q2 _ [p2 _ q 2 + Sq 2 


(4) 


fz |//p 2 . q2 + s 0 2 ) + 4p2q2 


8 1 


P \ \/\P 2 - q 2 + s o 2 ) 2 + 4p 2 q 2 - (p 2 - q 2 + s 0 2) + 


q\/l/(p 2 - q 2 + s 0 2j + 4p2q2 + p2 - q2 + Sq 2 


^2 ifp2 - q2 + s 0 2^ 2 + 4p2q2 


(5) 
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t 2 “ 


5 2 = 





J 

j/[l + So 2(p2 - <j2)j 

+ ^^2^2 + 1 + s 0 2 ( p 2 _ q 2 ) 


1 + s 0 2 

(p 2 - q 2 ) 

+ ^-s 0 ^p 2 q 2 

// 

1 + s 0 2 (p 2 - q 2 ) 

+ ^+ s cAp 2 q 2 - 

1 + s 0 2 (p 2 - q 2 ) 


( 6 ) 


(7) 


l/2 1 


1 + s o 2 (p 2 - <l 2 j 


+ 1 +s 0 ij -p 2 q 2 


Substitute t-^ and bp in the equations for A and a given below 
to determine Ap and and substitute T 2 and &2 to obtain A 2 

and 0 * 2 : 


a = 


1 + t 2 + 5 2 - ^(1 + t 2 + 6 2 ) - hr 2 - |l + (l - s Q ^j (' 




( 8 ) 


a = 


t 2 + 5 2 - A 


t 2 + \ _ 




(9) 


The elliptic functions E^]fo,s 0 2 J and E^l - p 2 m 2 j in the expres- 
sion for v (eq. (l)) may be obtained from tables (refs, 8 and 9) or by 
series expansion. 


Supersonic-Leading-Edge Wing 

The equation for the sidewash in the region below the right-hand 
panel of the wing and behind the Mach cone emanating from the wing apex 
(see appendix A) is 
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Y = 




rti / p 2 !!! 2 - 1 


1 + pm(p ^ 


tan 


-1 


. nk - ft |f - ft f 


tan' 


-1 


1 - pm(p £ 


- 11 - "I - ft I 


( 10 ) 


In the region behind the plane Mach wave off the wing leading edge 
and ahead of the Mach cone from the wing apex, the sidewash is given by 


v = 




1 /p^m 2 - 


(ID 


Sonic-Leading-Edge Wing 

When the wing leading edge is sonic, the sidewash expressions 
(eqs. (l) and (10)) simplify to 



p 
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RESULTS AND DISCUSSION 


Sidewash Charts 


By use of expressions derived herein, the sidewash for all points 
in the vicinity of infinite -chord triangular wings may be determined. 

In figures 2 to 7 the variation of the nondimens ional sidewash parameter 


with 


|3 ^ and 


H § for 


pm parameters of 0.4, 0.6, 0.8, 1.0, 1 . 25 , 


VooO, 

1«5> 1-75, 2.0, and 2.5 has been presented. Sidewash for pm values 
falling between those just enumerated may be obtained with sufficient 
accuracy in most cases by interpolation. In the succeeding discussion 

of the figures, the nondimens ional sidewash — — will be referred to 

Voo a 

simply as sidewash. 


Sidewash Variations With y and z 

By assu min g p and x to be fixed quantities in the sidewash vari- 
ations with p - and p it is possible to make a number of observa- 
x x 

tions on the variation of sidewash with vertical and spanwise distances. 
Another and more vivid method of picturing the flow is to construct con- 
tour plots of the sidewash, that is, lines of equal . — - . For this reason, 

V^a 

contour plots have been constructed for pm values of 0.8 and 1.5 and 
the results are shown in figures 8 and 9 * 

Sub sonic -leading -edge case . - On the surface of the subsonic-leading- 
edge wing (z = 0 ), the sidewash increases from a zero value at the wing 
plane of symmetry (y = 0) to infinity at the wing leading edge. (See 
fig. 2.) Outboard of the wing yet in the plane of the wing extended, 
the sidewash is zero. For a finite value of z, the sidewash increases 
with y from zero to a maximum, then decreases to zero at the Mach cone 
from the apex. For a given y value as z is increased (i.e., as z 
becomes more negative), figures 3 and 8 show that the magnitude of the 
sidewash decreases except in the region between the wing leading edge and 
Mach cone where there is an increase to a maximum from the zero value at 
the plane of the wing extended followed by a decrease. 

As the values of the parameter pm (for pm values less than 1) 
are decreased, the flow field below the wing remains similar in appearance 
and the magnitude of the sidewash decreases. 
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Sonic-leading-edge case .- When the wing leading edge is sonic, the 
sidewash is zero on the Mach cone everywhere except at the point where 
the wing leading edge and Mach cone coincide. It can be seen in fig- 
ure 4 that at this point the sidewash is infinite. 

Figure 5 presents the variation of sidewash with z for a number 
of y values and shows that for small values of y the sidewash decreases 
very slowly with z from its value at z = 0 over a large range of z 
values and drops rapidly near the Mach cone. For the larger values of y, 
the sidewash decreases more rapidly at the low negative values of z. At 
a fixed negative value of z, the sidewash increases from zero at the plane 
of symmetry to a maximum and then drops abruptly to zero at the Mach cone. 

Supersonic -leading-edge case .- In order to facilitate the discussion 
of the supersonic-leading-edge sidewash, a cross section of the wing and 
Mach cone system is given in the following sketch: 


i 



On the Mach cone emanating from the wing apex for the supersonic-leading- 
edge wing, the sidewash can be considered as having three different values. 
For that portion of the Mach cone that separates the region behind the Mach 
cone from the region between the Mach cone and the plane Mach wave off the 
leading edge, the sidewash has a constant value (eq. (11)). At the point 
where the plane wave off the leading edge is tangent to the Mach cone from 
the wing apex, the sidewash has half the constant value given by equa- 
tion (11). (See figs. 7(d) and 7(e).) Below the tangency point the side- 
wash has a zero value. 

The sidewash on the surface of the wing increases with y from a 
zero value at the wing plane of symmetry (y = 0) to a finite value at 
the Mach cone from the wing apex. For points on or below the wing out- 
side the Mach cone from the wing apex but within the plane Mach wave off 
the wing leading edge, the sidewash has a constant value (eq. (ll) ) . 
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The variation of sidewash with y (see figs. 6 and 9) for a given 
negative value of z, less than the negative z value of the tangency 
point, is very similar to the sidewash variation on the wing. For nega- 
tive values of z greater than that of the tangency point, the sidewash 
increases gradually with y until it reaches a maximum and then drops 
abruptly to zero at the Mach cone. 

Figures 7 and 9 show that the variation of sidewash with z for a 
fixed value of y is very slight over a large range of z values. Near 
the Mach cone, the sidewash drops abruptly to a zero value if the y value 
is less than that of the plane-wave tangency point and rises abruptly to 
a constant value if the y-coordinate is greater than that of the tangency 
point. 

Increasing pm (pm values greater than 1) causes a decrease in the 
maximum values of the sidewash. It is evident that as pm is increased 
the point of tangency of the leading-edge plane wave moves nearer to the 
wing symmetry plane and the region between the leading-edge plane wave 
and the apex Mach cone increases in size relative to the region behind 
the apex Mach cone. 


Sidewash Variations With x 

In order that the longitudinal or chordwise variations of sidewash 
may be conveniently plotted and more easily visualized, the semi-infinite 
triangular wing has been cut off at a distance x = c from the wing apex. 
The maximum distance behind the apex of the finite triangular wing that 
the sidewash can then be represented, by the equations derived herein, 
is 


x = c + pz 


or nondimens ionally 


x 

c 


1 + pm —J — 
b /2 


(13) 


where b/2 is the semispan of the finite-chord, triangular wing. Beyond 
the x value of equation ( 13 )> the flow field is affected by the finite- 
wing wake. 

The chordwise variation of the sidewash has been plotted in figures 10 
and 11 for two values of pm, a subsonic-leading-edge value of 0.8 and a 
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supersonic-leading-edge value of 1 . 5 $ three spanwise locations, = O. 25 , 

b/2 

0. 50 , and 0.75> and- three vertical heights, —7—= -0.1, -0.2, and -0.3* 

b/2 

For pm = 0.8, figure 10 shows that the sidewash rises rapidly to a maxi- 
mum from its zero value at the Mach cone and then decreases at a much more 
gradual rate to a finite value at the plane Mach wave off the trailing 
edge. As previously noted in the discussion of figures 2 and 3> the side- 
wash magnitudes increase in going from an inboard to an outboard station 
and decrease (though only slightly over the rear portion of the wing) as 
the vertical distance is increased. 

For the supersonic edge, pm = 1 . 5 , in figure 11, the most notable 
feature of the chordwise variations is the abrupt drop in the region 
immediately behind the Mach cone from the wing apex. In the rearward 
portions of the wing there is little change of sidewash with vertical 
height. This is consistent with previous observations. (See figs. 6, 

1 , and. 9.) 


CONCLUDING REMARKS 


Before an attempt can be made to calculate the aerodynamic forces 
acting on stores or missiles situated in the immediate vicinity of wings, 
it is necessary to have a detailed knowledge of the flow field near the 
wing. The present paper provides some of this needed information by 
presenting calculations of the lateral flow exterior to lifting, semi- 
infinite, triangular wings traveling at supersonic speeds. These solu- 
tions may also be used directly in determining the sidewash in some 
regions above and below finite wings. Charts are included from which 
the sidewash can be obtained, either directly or by interpolation, for 
all wing sweepback angles and supersonic Mach numbers. 

The methods used in deriving the sidewash are based on linearized 
theory; hence, the results are subject to its limitations. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., November 6, 1956. 
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APPENDIX A 

DETERMINATION OF SIDEWASH EXPRESSIONS 


The purpose of this appendix is to give the derivation of the side- 
wash equations. The derivations presented are carried out for p = 1 so 
that the sidewash and associated formulas given in subsequent sections 
are dependent only on m, y/x, and z/x. By use of the Prandtl-Glauert 
rule, which demonstrates that the flow-field velocities at any Mach number 

must depend on pm, p p and p the sidewash expressions may be 

rewritten to apply for all Mach numbers. (See eqs. (l), (10), (ll), and 
(12) in the main body of this paper.) 


Subsonic-Leading-Edge Case 

In reference 4, H. J. Stewart gives the following differential equa- 
tion (with a change in notation) for the complex sidewash velocity: 


dV 


- i 2 V oo aS 0 5 



(1 + £ 2 )d - 

U 2 + So 2 )(^So 2 



(Al) 


where 


s 


o 


1 - 


\i 


m 


m 


and 


i 


p + iq. 



(A2) 
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Integration of equation (Al) leads to 


V = v + iv 



The complex variables | and T) originate through variable changes 
necessary in the integration of equation (Al) and their relation to £ 
is given by 


6 


\| C 2 + So 2 


*1 



It is convenient to designate the real and imaginary parts of these new 
variables by 6-^, t 2 , and & 2 so that 


i = T 1 + i5 i 


and 


= t 2 - i5 2 

( T 5 , t 2 , and 6 2 in terms of y/x and z/x are given in the text 

by eqs. (2), (3), (4), (5), (6), and ( 7 )). The constant C in equa- 
tion (A3) is determined by the condition that the sidewash (real part of 
eq. (A3)) must be zero on the Mach cone. 
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Of interest here is the real part of equation (A3), that is, the side- 
wash. Thus, the problem remaining is to find the imaginary parts of the 
elliptic integrals E and F in equation (A3). With 

F (t ± iS, \j 1 - s 0 ^) = F r ± iFi 

and 

e(t ± 16 ,\/l - = Ep ± iEi 

the equation for the sidewash may be written 


v 



gVqjM Q 

s q 4 )e (\i 1 - m 2 ) 






+ C 


Reference 3 gives 


F i = F (\/a,s 0 2 ) 

Thus 

( P i) 1 * F (^o 2 ) 

K) 2 - P (^ s o 2 ) 


(Ak) 


(A5) 


(A6) 


and 


i6 
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In appendix B the imaginary part of the elliptic integral E with complex 
argument has been derived and yields 



[} - ( 1 " s o 4 ) A 1_ 

\i 


- 

ll - 

1 - 


■ s o 4 ) A l 

°1 


and 

( E i) 2 = " E (\[^2' s o 2 ) + F (>/W) + 


[l - (l - s 0 4 )A 2 _ 

\(cj^ ^1 - cr 2 \|l - s o 4 °2 

1 - 

H 

1 " s o^) A 2 

ct 2 


(A7) 


(a8) 


By substituting equations (A5)> (A6), (A7)> and (A8) into equation (A4), 
the sidewash is found to be 



p 
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In order to obtain the equation for v that is presented in the text, 
m (in the denominator of (eq. (A9)) has been replaced as follows by its 
equivalent in terms of s Q : 


m 


2s 


o 


1 


+ s 


2 

o 


Also, it has been found that the boundary conditions are satisfied with 
the constant of integration set equal to zero. 


Supersonic-Leading-Edge Case 

It has been demonstrated in reference 2 that the flow field exterior 
to a lifting or nonlifting semi-infinite triangular wing with supersonic 
leading edges is the sum of the flow fields of two lines of pressure 
sources coincident with the wing leading edges. Directly obtained in 
reference 2 by virtue of this result is the longitudinal perturbation 
velocity u for points in the flow field. The u velocity below the 
lifting wing and behind the Mach cone from the wing apex is given by the 
following equation: 


u 


cos- 1 * - . 

it \jm 2 - 1 \| (y - mx) 2 - z 2 (m 2 - l) 


COS 


X + 


EL 


7 


(y + mxr - ss ( 


m 


1) 


(A10) 


and for points between the Mach cone and the plane wave off the leading 
edge is given by 


-VqcP 1 
\] m 2 - 1 


(All) 


It should be noted that the positive square roots are to be taken in the 
denominator of the arc cosine terms of equation (A10) in conjunction with 
the principal values of the arc cosines. Equations for the vertical per- 
turbation velocity or downwash have been derived in reference 3. In the 
remainder of this section the lateral perturbation velocity, referred to 
herein as sidewash, will be obtained by the procedure used in reference 3. 
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By considering the relationship between the perturbation-velocity 
potential and the streamwise velocity component, 


dx 


= u 


expressions giving the perturbation-velocity potential 0 and lateral 
velocity component v may be written 


0 = J u(x' ,y,z)dx' 


( A12) 


= ^ = ^ J u(x',y,z)ax' 


(A13) 


The integration in equations (A12) and (A13) is along a line parallel to 
the x-axis from a point xq, where the potential 0 and velocity u are 
zero, to a point in the region where the velocity potential or sidewash 
is desired. 

Two possible paths of integration are shown in figure 12 for points 
behind the apex Mach cone. One path intersects both the plane wave and 
Mach cone whereas the other intersects only the Mach cone. Mathematically, 
it can be shown that the first of these integration paths corresponds to 
the condition that x < my and the second, to the condition that x > my. 
Consider the first of the two paths of integration. The velocity potential 
at a point x from equations (A10), (All), and (A12) is given as follows: 


a , - v »“ n 

n \ly 2 +z 2 

/ , dx' + i 

f x 

/ cos”^- - 

x' - my 

r~H 

J y z\|m 2 -l 
m m 


j{ y - mx' ) 2 - z 2 (m 2 - 1) 


r 


x ! + my 

\|(y + mx’ )2 - z 2 (m^ - l) 


dx ' 


1 


(Al4) 
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Performing the first two integrations and then differentiating with respect 
to y gives rise to 


V^om 

1 

1 

>> 

y 

n 


\Jy2 + m\fm^ - 1 

Jy 2 + z 2 






mx' 2 - x'y - mz 2 


\|y2 +z 2 \| x '2 _ y 2 _ z 2 j^y _ mxi )2 _ z 2( m 2 _ 


dx ' 


\ 



The first two terms in equation (A15) are obtained from the first integral 
of equation ( Al4 ) . The third and fourth terms result from differentiating 
the first arc cosine integral in equation (Al4) and the last term in equa- 
tion (A15) comes from differentiating the second arc cosine term of equa- 
tion (Al4). Note that the first and third terms of equation (A15) cancel. 
The first integral in equation (A15) when evaluated (remember when sub- 
stituting the lower limits that x < my) gives 


1 1 

tan" 1 x^JSSL 

+ \/ni^ - 1 COSh"^ ^ ^ 1 

r— 1 

f 

OJ 

£ 

( \|m 2 - 1 \|x 2 - y 2 - z 2 

l 2 + ^ 2 


(Al6) 

The evaluation of the second integral in equation (A15) may be obtained 
from the evaluation of the first integral in equation (A15) (expres- 
sion (Al6) ) by replacing y by -y and changing the sign before the 

■4 term. This procedure yields 


mum 


tan 


-1 


x + my 


(jin 2 " 


- 1 


f2 „2 


\fm 2_ 


- 1 cosh 


-1 


x - y - z 




x 

2 2 
y + 



(A17) 
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When expressions (Al6) and (A17) are substituted for the integrals in 
equation (A15), the sidewash behind the apex Mach cone for the case where 
x < my is seen to be 


v = -- 


Vco* 


tan 


-1 x - my 


- tan 


-1 x + my 


it \fi 


m 2 - i 


\/m 2 - 1 \[: 


2 , r 2 „2 


tot - 1 yx - y - z 


\jm 2 - 1 ]fx^ - y 2 - z 2 

(A18) 


When the path of integration to the point x intersects only the 
Mach cone, the sidewash from equations (A10) and (A13) may be written as 


v = - 


3 

/ am ^ 

00 O 

f X -1 

/ COS 

x 1 - my 

M 

m2 - 1 ^ ' 

\jy 2 +z 2 

J(y - mx^ 2 - z 2 (m 2 - l) 


dx' + 


cos 


-1 


i y 2 +z 2 


x.^ + my 


\[{y + mx x ) 2 - z 2 (m 2 - l) 


dx' 


(A19) 


Carrying out the differentiation in equation (A19) yields 


V 


Vjxm 

it 


W 


2 2 
y +z 


mx 


2 


x'y - mz 


2 



|(y - mx’ ) 2 - 




dx’ 





y 
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which when integrated with x > my results in the following equation: 


V„a 




n \|m 2 - 1 


tan 


-1 


x-ray 


y \|m 2 - 1 \jx 2 - y 2 - z 2 


+ \)ra 2 - 1 cosh 1 


\Jy 2 + z 2 


tan 


-1 


x + ray 

l/m 2 - 1 \/x 2 - y 2 - z 2 




1 cosh" 





tan 1 



x -f ray 

1 |/x 2 - y 2 - z 2 


( A20) 

It can be seen from equations (Al8) and (A20) that the sidewash behind 
the Mach cone is given for all points by the same expression. 

The sidewash in the region between the plane Mach wave off the leading 
edge and the Mach cone below the right wing panel is given by 


v 


(^1 ay 



dx' 


m m 


which, when the indicated operations are performed, becomes 
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Sonic-Leading-Edge Case 

The sonic -leading -edge sidewash is most simply obtained from equa- 
tion (A20). Setting m equal to 1 in equation (A20) gives rise to an 
indeterminancy which may easily be evaluated to yield the required sidewash 


v 



(A22) 
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APPENDIX B 

EVALUATION OF IMAGINARY PART OF AN INCOMPLETE ELLIPTIC 
INTEGRAL OF THE SECOND KIND WITH A COMPLEX ARGUMENT 


In reference 3 the real part of the incomplete elliptic integral 
E(T+i5,k) and the real and imaginary parts of the elliptic integral 
F(T+i5,k) have been derived. The purpose of this appendix is to obtain 
the imaginary part of the elliptic integral E(T±i5,k) which is necessary 
in determining the sidewash. (See eq. (A4).) 

The incomplete elliptic integral of the second kind with a complex 
argument may be written as 


where 


E(T±i6,k) = f 

J r 


T±i6 


/l - k 2 | 2 

\/l - I 2 


d| 


k = 


v/l 


- s 


k 

o 


(Bl) 


When the Jacobian transformation 


T ± 18 = i = sn(|i,k) = sn(n r ±ni,k) 


(B2) 


is introduced, equation (Bl) becomes 


E( T ±i5,k) 


-L 


F +iF- 
r u i 


Q 

dn^(|i,k)d)! 


(B3) 


Evaluation of the real and imaginary parts of the upper limit of equa- 
tion (B3), that is, 


F r ± iF^ 





is given in reference 3. (See also eqs. (A5) and (A6) . ) 
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As pointed out in reference 3, the integration in equation (B3) is 
most conveniently performed in two steps: first, along the real axis 

to F r , and then parallel to the imaginary axis from F r to F r ± iFj_. 

Thus, equation (B3) becomes 


E(r±iS,k) 



dn 2 (p r ,k)dp r 




i 



dn 


(F r ±ip i ,k)dp.i 


(BU) 


The first integral in equation (b 4) is real and may be excluded from 
further consideration. It is convenient in determining the imaginary 
part of the second term to expand the delta-amplitude elliptic function 
with the aid of the following formula given in reference 8, page 24: 


dn(F r ,k)cn(pi,k' )dn(pi,k' ) _ ik 2 sn(F r ,k)cn(F r ,k)sn(pi,k' ) 

dn(F r ±ipi,k) = p 5 + p“ 0 “ 

1 - sn (pi,k')dn (F r ,k) 1 - sn (pi,k')<hi (F r ,k) 

With the use of the relationship in the second integral of equation (B4), 
the imaginary part of E(T±i5,k) is given by 

dn 2 (F r ,k)cn 2 (pi,k')dn 2 (pi,k' ) 

[l - sn 2 (pi,k' )dn 2 (F r ,k)] 2 

k^sn^Fj^k^n^Frjk^n^iaj^k' ) 

— r f<3Pi (B3 

[l - sn 2 (pi,k' )dn 2 (F r ,k)J 2 


r? i 

I.P.[E(T±i5,k)] = ± / 

J f) 


By use of the integration formulas given on pages 218 and 219 of refer- 
ence 8, the imaginary part of E(T±i5,k) in terms of Jacobian elliptic 
functions becomes 


I.P. [ E(r±iB,k)J 


+ 


F^sn Fi,k') 


E( sn Fi,k') + 


dn 2 (F r ,k)sn(Fi,k' )cn(Fj.,k' )dn(Fj_,k' ) 


1 - dn 2 (F r ,k)sn 2 (Fi,k' ) 


(B6) 
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From the change in variable of equation (B2) , that is, 

t ± iB = sn(F r ±Fi,k) 


and by defining A and a as follows: 


A = sn 2 (F r ,k) 
a = sn 2 (Fi,k' ) 


(B7) 

(B8) 


reference 3 has shown that t and 6 may he related to A and a by 


t 2 _ A(l - k ,2 cr) 


(l - a + k 2 Acr) 2 


(B9) 


2 a(l - a) (l - A) (l - k 2 A) 
o = 

(l - a + k 2 Acr) 2 


(BIO) 


and 


[l + t 2 + 6 2 - + t 2 + 5 2 ) 2 - 4t 2 ] {l + k 2 (-r 2 + 6 2 ) - J[l + k 2 (r 2 + 6 2 ) 2 - 4k 2 T 2 } 


hk 2 T 2 


(Bll) 


t 2 + 5 2 - A 

t 2 + 8 2 - A - [Ak 2 (r 2 + 5 2 ) - l]] 


(B12) 


Other useful relationships which may be formed by use of equations (B7) 
and (B8) are 


cn(F r ,k) = ^/l - A 


(B13) 


dn(F r ,k) = 


k 2 A 


(B14) 
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cn(Fi,k') = <Jl - a 
dn(Fi,k') = \fl - k'^a 


(B15) 

(Bl6) 


By use of the relationships given in equations (B7), (b 8 ), (B13), (Bl4), 
(B15), and (Bl 6 ), the imaginary part of the elliptic function E(T±iS,k), 
defined by equation (b 6 ), can now be written in terms of A, a, and k 
as 


I.P.[E(r±i5,k)] = ± F^k’j - E k' 

(l - k 2 A) s/g \/l~- a \/l - k l2 o 
1 - (l - k 2 A)a 


In equations (A 6 ) and (A 7 ), k and k' have been replaced by their 


equivalents 




s 0 4 and s 0 £ 
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Figure 1.- Sketch of supersonic-leading-edge wing showing system of axes used in analysis and 
the positive directions of the flow velocities and angle of attack. 
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(a) pm = 0.4. 

Figure 2.- Variation of the subsonic-leading-edge sidewash with pi for a range of p £ 

and Pm values. 
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Figure 2.- Continued. 



Figure 2.- Concluded 
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(a) (3m = 0.4. 

Figure 3*- Variation of the subsonic-leading-edge sidewash with 0— for a range of (3*4 

X x 

and pm values . 
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Figure 3«- Continued. 
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(c) Pm = 0.8. 


Figure 3«- Concluded. 
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Figure k.- Variation of the sonic-leading~edge (|3m = 1) sidewash with (3^ 

for a range of (3^ values. 
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Figure 5 .- Variation of the sonic -leading-edge (pm = 1) sidewash with p^- for a range of 

P^ values. 
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(a) pm = 1 . 25 . 

Figure 6.- Variation of the conical part of the supersonic -leading-edge 
sidewash with p| for a range of p| and pm values. 
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(b) pm = 1.5. 
Figure 6.- Continued. 
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Figure 6.- Continued 
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(d) pm = 2.0. 
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(e) pm = 2.5. 


Figure 6.- Concluded. 
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(a) pm = 1 . 25 . 

Figure 7«- Variation of the conical part of the supersonic-leading-edge 
sidewash with p£ for a range of p^ and pm values. 
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Figure 7»- Continued. 
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(c) pm = I.75. 
Figure 7.- Continued. 
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(a) pm = 2 . 0 . 



(e) pm = 2 . 5 . 


Figure 7 .- Concluded. 
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Figure 8.- Sidewash contours for a subsonic-leading-edge lifting tri- 
angular wing of infinite chord. ( 3 m = 0.8. 
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Figure 9.- Sidewash contours for a supersonic-leading-edge lifting triangular wing of infinite 

chord, pm = 1.5 • 
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(a) JL. = 0 . 25 . 
b/2 
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Figure 10.- Chordwise variation of sidewash at three spanwise positions 
and three vertical heights for a subsonic leading edge. Pm = 0.8. 
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(a) -2- = 0.25. 

b/2 



(b) = 0.50. 

b/2 



Figure 11.- Chordwise variation of the conical part of the sidewash at 
three spanwise positions and three vertical heights for a supersonic 
leading edge. Pm = 1.5* 


NACA - Langley Field, Va. 



Figure 12.- Two paths of integration for obtaining the velocity potential and associate notation. 
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